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Masses and widths of the axial- vector charmed mesons Di(2420), Di(2430), D s i(2536), and 
D s i(2460) are calculated nonperturbatively in the Resonance-Spectrum- Expansion model, by cou- 
pling various open and closed meson-meson channels to the bare J p = 1 + cq (q — u, d) and cs 
states. The coupling to two-meson channels dynamically mixes and lifts the mass degeneracy of 
the spectroscopic 3 Pi and 1 Pi states, as an alternative to the usual spin-orbit splitting. Of the two 
resulting S'-matrix poles in either case, one stays very close to the energy of the bare state, as a 
quasi-bound state in the continuum, whereas the other shifts considerably. This is in agreement 
with the experimental observation that the Di(2420) and D s i(2536) have much smaller widths than 
one would naively expect. The whole pattern of masses and widths of the axial- vector charmed 
mesons can thus be quite well reproduced with only two free parameters, one of which being al- 
ready strongly constrained by previous model calculations. Finally, predictions for pole positions of 
radially excited axial-vector charmed mesons are presented. 

PACS numbers: 14.40. Lb, 13.25.Ft, 11.80.Gw, 11.55.Ds, 12.40.Yx, 



I. INTRODUCTION 

The axial- vector (AV) charmed mesons Di(2420) and 
D sl (2536) pQ have the puzzling feature that their de- 
cay widths are much smaller than one would expect 
on the basis of their principal S'-wave decay modes. 
Namely, the Z?i(2420) decays to D*ir (possibly also in 
a D wave), with a phase space of more than 270 MeV, 
but has a total width of only 20-25 MeV pQ. On the 
other hand, the D s i(2536) decays to D* K in S and D 
wave with a phase space of about 30 MeV, resulting in 
an unknown tiny width <2.3 MeV, limited by the ex- 
perimental resolution pQ. The discovery of the miss- 
ing two AV charmed mesons, namely the very narrow 
D s i(2460) and the very broad L>i(2430), first observed 
by CLEO [2] and Belle [3J, respectively, completed an 
even more confusing picture. While the tiny width of 
the D s i(2460) can be easily understood, since this meson 
lies underneath its lowest Okubo-Zweig-Iizuka-allowed 
(OZIA) and isospin-conserving decay threshold, the huge 
Di(2430) width, in D*ir, is in sharp contrast with that of 
the L>i(2420). Moreover, the L> s i(2460) lies 76 MeV be- 
low the L> si (2536), whereas the £>i(2420) and L>!(2430) 
are almost degenerate in mass, if one takes the central 
value of the latter resonance. 

Quark potential models, with standard spin-orbit split- 
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tings, fail dramatically in reproducing this pattern of 
masses. For instance, in the relativized quark model [3] 
the cs state that is mainly 3 Pi comes out at 2.57 GeV, as- 
suming the already then well-established Z? sl (2536) to be 
mostly 1 P\ , though with a very large mixing between 3 P\ 
and 1 P\ . Reference [1] similarly predicted a too high mass 
for the dominantly 3 P\ state in the cq (q = u, d) sector, 
viz. 2.49 GeV. In the chiral quark model for heavy-light 
systems of Ref. [5] , the result for the mainly 3 Pi cq state is 
also 2.49 GeV, while the discrepancy is even worse in the 
cs sector, with a prediction of 2.605 GeV for the mostly 
3 Pi state, now with a small mixing in both sectors. 

More recently and after the discovery of the D sl (2460) 
(and Di(2430)), chiral Lagrangians for heavy- light sys- 
tems (see e.g. Refs. [5HS]) have been employed in order 
to understand the masses of the AV charmed mesons, in 
particular the mass splittings with respect to the vector 
(V) mesons with charm D* and D* , respectively. Ref- 
erence [7] analyzed in detail the curious experimental [l! 
observation that the AV-V mass difference is considerably 
larger in the charm-nonstrange sector than in the charm- 
strange one, which is not predicted by typical quark po- 
tential models [3J|S]. The same discrepancy applies to the 
scalar-pseudoscalar mass difference in either sector [TJ [7] . 
In Ref. [7] , the problem was tackled by calculating chiral 
loop corrections, but the result turned out to be exactly 
the opposite of what is needed to remove or alleviate the 
discrepancy. 

An alternative approach to the AV charmed mesons 
is by trying to generate them as dynamical resonances 
in chiral unitary theory [10] . Indeed, in the latter pa- 
per, describing AV mesons in other flavor sectors as 
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well, several charmed resonances were predicted, includ- 
ing the £»i(2420), A. (2430), D sl (2536), and L» s i(2460), 
with reasonable results, though the cq states came out 
about 100 MeV off. However, dynamical generation of 
mesonic resonances, including the ones that are com- 
monly thought to be of a normal quark-antiquark type, 
may give rise to interpretational difficulties, besides pre- 
dicting several genuinely exotic and so far unobserved 
states [TU]. Dynamically generated AV charmed as well 
as bottom mesons can be found in Ref. [11 , too. 

Finally, in Ref. [Hj a coupled-channel calculation of 
positive-parity cs and bs was carried out in a chiral quark 
model, similar to our approach in its philosophy, and 
with results for the Z? s i(2536) and D s i(2460) close to 
the present ones (also see below). 

II. RESONANCE-SPECTRUM EXPANSION 

In the present paper, we employ the Resonance- 
Spectrum Expansion (RSE) to describe the AV charmed 
mesons. The RSE model has been developed for meson- 
meson (MM) scattering in non-exotic channels, whereby 
the intermediate state is described via an infinite tower of 
s-channel qq states 13J . For the spectrum of the latter, 
in principle any confinement potential can be employed, 
but in practical applications, a harmonic oscillator (HO) 
with constant frequency has been used, with excellent 
results 13J . Recent applications of the RSE concern the 
(£(2170) Q3] and A(3872) [15] resonances. 

In order to account for the two possible spectroscopic 
channels 3 Pi and 1 P\ contributing to a J p = 1 + state 
with undefined C-parity, we couple both qq channels to 
the most important meson-meson channels. The result- 
ing fully off-energy-shell RSE T matrix reads [13TH5] 

T^'^ipup'^E) = -2A 2 7MOT4(Pir )x 

N 

£ K im { [i - nn]- 1 }^ 4 (p» yfe> , (i) 

m— 1 

with the diagonal loop function 

ttij(kj) = -2iX 2 fi j k j r j J L .(k j r )h^j :i (k j ro)Si j , (2) 
and the RSE propagator 

S=0,ln=0 & fi « 

Here, A is an overall coupling, ro is the average distance 
for decay via 3 Po quark-pair creation, is the dis- 

crete energy of the n-th recurrence in the qq channel 
with spin S, g^ s n ^ is the corresponding coupling to the 
j-th MM channel, ^ the reduced mass for this channel, 
Pi the off-shell relative momentum, Li the orbital an- 
gular momentum, and j l L , and (kjro) the spherical 



Bessel and Hankel functions of the first kind, respectively. 
Note that /Zj, Pi, and the on-energy-shell relative mo- 
mentum ki are defined relativistically. Also notice that 
the infinite sum over the higher recurrences converges 
very fast, so that it can be truncated after 20 terms in 
practical calculations. The S matrix is finally given by 
S% uLs) {E) = 1 + 2iT^ i ' L '\k i) k j ;E). 



III. OZI-ALLOWED CHANNELS FOR AV 
CHARMED MESONS 

Now we describe the physical AV charmed resonances 
by coupling bare 3 Pi and 1 P\ cq, cs channels to all 
OZI-allowed ground-state pseudoscalar-vector (PV) and 
vector- vector (W) channels. It is true that there are also 
relevant pseudoscalar-scalar (PS) channels (in P-wave), 
most notably £>/ (600) and D5(2400)tt Q] in the AV cq 
case, and DKq (800) for cs. These will contribute to the 
observed pQ Dttit and DttK decay modes, respectively. 
Now, we have recently developed |15j an algebraic proce- 
dure to deal with resonances in asymptotic states whilst 
preserving unitarity. However, the huge widths of the 
D2(2400), / (600), and if *(800) resonances may lead to 
fine sensitivities that will tend to obscure the point we 
want to make, apart from the fact that there will also 
be nonresonant contributions to the Drnr and DnK final 



TABLE I: Included meson-meson channels for Di(2420) and 
Di(2430), with ground-state couplings squared [TS], orbital 
angular momenta, and thresholds in MeV. For r\ and jj' , a 
pseudoscalar mixing angle of 37.3° 114^ is used. 



Channel 


(fl(S=l,n=0)) 


(5(S=0,n=0)) 


2 L 


Threshc 


D*TT 


0.02778 


0.01389 





2146 


D*n 


0.03472 


0.06944 


2 


2146 


D*rj 


0.00524 


0.00262 





2556 


D*r] 


0.00655 


0.01310 


2 


2556 


D* 3 K 


0.01852 


0.00926 





2608 


D*K 


0.02315 


0.04630 


2 


2608 


Dp 


0.02778 


0.01389 





2643 


Dp 


0.03472 


0.06944 


2 


2643 


Du 


0.00926 


0.00463 





2650 


Duj 


0.01157 


0.02315 


2 


2650 


D*p 





0.01389 





2784 


D*p 


0.01042 


0.06944 


2 


2784 


D'oj 





0.00463 





2791 


D*uj 


0.03472 


0.02315 


2 


2791 


D S K* 


0.01852 


0.00926 





2862 


D S K* 


0.02315 


0.04630 


2 


2862 


D*r) 


0.00402 


0.00201 





2996 


D*vl 


0.00502 


0.01004 


2 


2996 


D*K* 





0.00926 





3006 


D*K* 


0.06944 


0.04630 


2 


3006 



TABLE II: As Table[E) but now for D 3 i(2536) and D s i(2460). 



Channel 


(§(S=l,n=0)) 


G?(S=0,n=0)) 


* L Threshc 


D*K 


0.03704 


0.01852 





2504 


D*K 


0.04630 


0.09259 


2 


2504 




0.00803 


0.00402 





2660 


Din 


0.01004 


0.02009 


2 


2660 


DK* 


0.03704 


0.01852 





2761 


DK* 


0.04630 


0.09259 


2 


2761 


D*K* 





0.01852 





2902 


D*K* 


0.01389 


0.09259 


2 


2902 


D a <j> 


0.01852 


0.00926 





2988 


D 3 4> 


0.02315 


0.04630 


2 


2988 


DU 


0.01048 


0.00524 





3069 


DW 


0.01310 


0.02621 


2 


3069 


D*A 





0.00926 





3132 


D* s <j> 


0.06944 


0.04630 


2 


3132 



states. So we restrict ourselves to the open and closed 
PV and VV channels in the present investigation, but we 
shall further discuss this issue below. The here included 
channels for cq and cs are given in Tables [T| and |H} respec- 
tively, together with the corresponding orbital angular 
momenta, threshold energies, and ground-state couplings 
squared {g^ s=1 ^ n=0 ) 2 , where S = 1(0) refers to the 
3 Pi ^-Pi) quark- antiquark component. In Appendix [Aj 
we show in more detail how the ground-state coupling 
constants in Tables [I] and [TT] depend on the isospin and 
J quantum numbers of the various meson-meson chan- 
nels. The latter squared couplings, computed in the very 
general framework of Ref. [16] , must be multiplied by 
(n + l)/4" for L = and by (2n/5 + l)/4 n for L = 2, so 
as to obtain the couplings for the radial recurrences n in 
the RSE sum of Eq. ([3} . Note that the scheme of Ref. Ql] 
employs overlaps of HO wave functions for the original 
qq pair, the 3 P qq pair created out of the vaccuum, and 
the outgoing mesons. This allows to rigorously calculate 
the coupling constants of all excited states as well, in 
contrast with approaches using combinations of Clebsch- 
Gordan coefficients only. Nevertheless, our ground-state 
couplings are identical to the usual ones in practically 
all situations, including the present one. Finally, a sub- 
threshold suppression of closed channels is used just as 
in Ref. HI]. 

The energies of the bare AV cq and cs states we deter- 
mine, as in previous work (see e.g. Refs. |141 115) ). from 
an HO spectrum. The corresponding constant oscillator 
frequency and the constituent masses of the charmed, 
strange, and nonstrange quarks are also kept completely 
unchanged at the values lj = 190 MeV, m c = 1562 MeV, 
m s = 508 MeV, and m n = 406 MeV [12 HI]. This yields 
masses of 2443 MeV and 2545 MeV for the bare AV cq 
and cs states, respectively, which are very close to values 
found in typical single-channel quark models [JJ [5] . 
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FIG. 1: Di(2430) pole trajectories as a function of A, for 
ro = 3.2-3.5 GeV -1 (left to right). Solid curve and bullets 
correspond to ro = 3.40 GeV -1 and A = 1.30, respectively. 



IV. QUASI-BOUND STATES IN THE 
CONTINUUM AND OTHER POLES 

Next we search for poles in the S matrix. Start- 
ing with the cq case, we choose r in the range 3.2- 
3.5 GeV -1 (0.64-0.70 fm), which is in between the values 
of 2.0 GeV" 1 [T5] for an AV cc system and 4.0 GeV" 1 [H] 
for vector ss states. In Fig.[TJ we plot several pole trajec- 
tories in the complex E plane as a function of the overall 
coupling A. We see that this pole rapidly acquires a large 
imaginary part, whereas the real part changes consider- 
ably less, especially in the range r = 3.3-3.5 GcV -1 , 
making it a good candidate for the broad Di(2430) res- 
onance. For A = 1.30 and r — 3.40 GeV -1 , the pole 
comes out at (2430 — i x 191) MeV, being thus fine-tuned 
to the experimental mass and width pQ. However, there 
should be another pole in the S matrix, since there are 2 
quark-antiquark channels and more than 2 MM channels. 
From the structure of the T matrix in Eqs. ([l]{3]), one can 
algebraically show that the number of poles for each bare 
state is equal to min(N q q, Nmm), besides possible poles 
of a purely dynamical nature. Indeed, another pole orig- 
inating from the bare cq state is encountered, with its 
trajectories depicted in Fig. [2] Quite remarkably, this 
pole moves very little, acquiring an imaginary part that 
is a factor 55 smaller than in the £>i(2430) case, for the 
values A = 1.30 and r — 3.40 GeV -1 (see solid lines and 
bullets in both figures). So this resonance, with a pole 
position of (2439 - ix 3.5) MeV, almost decouples from 
the only open OZIA MM channel [T7], viz. D*n, repre- 
senting a quasi-bound state in the continuum (QBSC) 



ReE (GeV) 
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FIG. 2: Di(2420) pole trajectories as a function of A, for 
ro = 3.3-3.5 GeV -1 (left to right). Solid curve and bullets 
correspond to ro = 3.40 GeV -1 and A = 1.30, respectively. 



[IT]. Moreover, it is a good candidate for the i?i(2420), 
though its width of roughly 7 MeV is somewhat too small 
and its mass 16 MeV too high. These minor discrepan- 
cies may be due to the neglect of the PS channels, with 
broad resonances in the final states, as suggested above. 

Nevertheless, these encouraging results might be partly 
due to a fortuitous choice of the parameters A and r . 
Therefore, we now check the cs system, thereby scal- 
ing ro and A with the square root of the reduced quark 
mass (see Ref. [17] , Eq. (13)), so as to respect flavor in- 
dependence of our equations, which yields the cs val- 
ues r = 3.12 GeV -1 and A = 1.19. The ensuing cs 
pole trajectories are depicted in Fig. [3j but now for 
ro = 3.12 GeV -1 only. Thus, for A = 1.19, the strongly 
coupling state comes out at 2452 MeV, i.e., only 7.5 MeV 
below the D sl (2460) mass, with a vanishing width, as the 
pole ends up below the lowest OZIA channel. As for the 
cs QBSC, it indeed shifts very little from the bare state, 
settling at (2540 - i x 0.7) MeV, i.e., only 5 MeV above 
the -Dai (2536) mass, and having a width fully compatible 
with experiment pQ. 

Besides the above ground-state AV charmed mesons, 
the present model of course also predicts higher recur- 
rences of these resonances. However, due caution is nec- 
essary so as to account for the most relevant open and 
closed decay channels at the relevant energy scales. Now, 
the first radially excited HO levels of the 3 Pi/ 1 Pi cn and 
cs states lie at 2823 MeV and 2925 MeV, respectively, 
which allows the corresponding resonances to be rea- 
sonably described by the channels included in Tables [I] 
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[IlJ Thus, we find again 4 poles, tabulated in Table III, 
together with those of the ground-state AV charmed 

TABLE III: Poles of ground-state (rt = 0) and first radially- 
excited (n = 1) AV charmed mesons. Parameters: A = 1.30 
(1.19) and r = 3.40 (3.12) GeV -1 , for cq (cs) states. 

Quark Content Radial Excitation Pole in MeV 



cq 





2439 - i X 3.5 


cq 





2430 - i x 191 


cs 





2540 - i X 0.7 


cs 





2452 - ixO.O 


cq 


1 


2814 -i x 7.8 


cq 


1 


2754 - i x 47.2 


cs 


1 


2915 - i X 6.7 


cs 


1 


2862 - i x 25.7 



mesons. For the radially excited states, we observe a sim- 
ilar pattern as for the ground states, namely two poles 
that remain close to the bare HO levels, whereas two 
other poles shift considerably. Note, however, that the 
difference is not as dramatic as in the n — case. This 
may be due to the fact that several decay channels are 
open now. As for a possible observation of the here pre- 
dicted 2 Pi states, no experimental candidates have been 
reported so far. Namely, in the nearby cq mass region, 
the two listed Q] resonances L>(2600) and D(2750) [I] 
both decay to D*ir and Dn, which excludes an AV as- 
signment. 
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FIG. 3: D s i(2460) (dashed) and D s i(2536) (solid) pole tra- 
jectories as a function of A, for ro = 3.12 GeV -1 . Bullets 
correspond to A = 1.19; vertical line shows D* K threshold. 
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Concerning the cs sector, the only listed pQ state 
around 2.8-2.9 GeV is the L>*j(2860) [E], with natu- 
ral parity and so not an AV, decaying to D*K and DK, 
which makes it a good candidate for the 2 3 p2 state, pos- 
sibly overlapped by the 2 3 Po [19 . Note that the lower of 
our two predicted 2 Pi resonances also practically coin- 
cides with the D* j (2860), both in mass and width. This 
may be a further indication that the D*j(2860) structure 
corresponds to more than one resonance only. 

To conclude this section, we study — for the cq sys- 
tem — the dependence of the lowest-lying poles on the 
number of included quark-antiquark and MM channels. 
In Tabic IWJ besides the £>i(2420) and L>i(2430) poles 



TABLE IV: Poles of AV cq mesons, for different sets of in- 
cluded channels. Parameters: A = 1.30, ro = 3.40 GeV -1 . 

cq channels MM channels Pole 1 (MeV) Pole 2 (MeV) 





20 


2430- 


i x 191 2439 


-i x 3 




2 


2402 - 


i x 36 2441 


— i x 1 


S P1+ X P1 


1 


2431 - 


i x 39 




3 A 


20 


2409- 


i x 65 




"Pi 


20 


2425 - 


i x 96 





resulting from the full calculation, with the 20 MM chan- 
nels from Table |TJ we first give the pole positions for the 
cases that only 2 (D*tt, L = 0, 2) or 1 (D*n, L = 0) 
MM channels are included. The last two poles then cor- 
respond to calculations with the full 20 MM channels but 
only one quark-antiquark channel, viz. 3 Pi or 1 P\. Notice 
that only one pole is found when the number of quark- 
antiquark or MM channels is equal to 1. This confirms 
our above conjecture that the number of poles for each 



bare HO level is given by m.m(N q q, 



MM 



the decay couplings in Tables [T] and [TTJ dispensing with 
the usual S ■ L splitting. Also note that the occurrence 
of (approximate) bound states in the continuum for AV 
charmed mesons had already been conjectured by two of 
us |17j . based on more general arguments. 

The puzzling discrepancy between the AV-V mass 
splittings in the cq and cs sectors is resolved in our cal- 
culation by dynamical, nonperturbative coupled-channel 
effects. A similar phenomenon we have observed before 
[21] for the (2300-2400) [T] resonance, and may be re- 
lated to an effective Adler-type zero [22] in the D*tt and 
Dir channels in the AV and scalar cn cases, respectively, 
owing to the small pion mass. 

Summarizing, we have reproduced the whole pattern of 
masses and widths of the AV charmed mesons dynami- 
cally, by coupling the most important open and closed 
two-meson channels to bare cq and cs states contain- 
ing both 3 Pi and 1 P\ components. The dynamics of the 
coupled-channel equations straightforwardly leads to one 
pair of strongly shifted states and another pair of QBSCs. 
Ironically, the state that shifts most in mass, namely the 
D s i(2460), ends up as the narrowest resonance. This em- 
phasizes the necessity [53] to deal with unquenched me- 
son spectroscopy in a fully nonperturbative framework. 

One might argue that these conclusions will depend on 
the specific model employed. Admittedly, our numerical 
results could change somewhat if slightly different bare 
masses for the AV charmed mesons were chosen, non-iS 1 - 
wave decay channels were included as well, or a different 
scheme was used to calculate the decreasing couplings of 
the higher recurrences. Nevertheless, we are convinced 
the bulk of our results will not change, most notably 
the appearance of QBSCs and the large shifts of their 
partner states, as the almost inevitable consequence of 
exact nonperturbative coupled-channel dynamics. 
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In the foregoing, we have managed to rather accu- 
rately reproduce the masses and widths of the D\ (2420) , 
L>i(2430), P s i(2536), and D s i(2460) with only 2 free pa- 
rameters, one of which is already constrained by previ- 
ous model calculations, as well as by reasonable estimates 
for the size of these mesons. Crucial is the approximate 
decoupling from the continuum of one combination of 
3 Pi and l P\ components, which amounts to a mixing 
angle close to 35°. Namely, if we express a QBSC as 
IQBSC) = -sin(9| 3 Pi) +cos0| 1 Pi), it decouples from 
the L — D*tt channel (for cq) or D*K channel (for cs), 
if 6 = arccos y/2/3 « 35.26° (see Tables |lJH. Inclusion 
of the other, practically all closed, channels apparently 
changes the picture only slightly in our formalism. This 
result is in full agreement with the findings in Ref. [T2"] . 
However, in the present approach this particular mixing 
[20j comes out as a completely dynamical result, and is 
not chosen by us beforehand. Moreover, the bare-mass 
degeneracy of 3 Pi and X P\ states is adequately lifted via 
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Appendix A: Three-meson couplings 

The ground-state couplings in Tables [T] and |TT] are ob- 
tained by multiplying the isospin recoup ling s given in Ta- 
ble with the J PC couplings in Table fvi] for an OZIA 
process Ma — > Mb + Mc based on 3 P qq creation \Tjg\ . 
For clarity, we represent here all couplings by rational 
numbers. Note that r\ n and r) s in Tabl e [V] stand for the 
pseudoscalar 1 = states (uu + dd)/\ / 2 and ss, respec- 
tively. Then, we get the couplings to the physical r\ and 
r/ mesons by applying a mixing angle — in the flavor 
basis — of 41.2°, as in Ref. [21], 2nd paper. For the ui 
and <p we assume ideal mixing. 
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TABLE V: Squared isospin recouplings for the 3-meson pro- 
cess Ma — > Mb + Mc , with Ma = cs or cq. 



M A 


Mb 


Mo 


o! 


D sl 


Ds,D* a 


Vs,<t> 


1/3 


D sl 


D, D* 


K,K* 


2/3 


D 1 


Ds,d; 


K, K* 


1/3 


D 1 


D,D* 


n,p 


1/2 


Di 


D,D* 


fin, tO 


1/6 



TABLE VI: Squared ground-state coupling constants for the 
3-meson process M A -> M B + Mc, with J PC (M A ) = l +± , 
and Ma, Mb belonging to the lowest pseudoscalar or vector 
nonet. 



PC (M A ) 


J PC (M B ) 


J PC (M C ) L 


M B M C 


Sm b 


2 

Mc 5(n=0) 


1++ 


o-+ 


1~~ 





1 


1/18 


1++ 


o-+ 


1~~ 


2 


1 


5/72 


1++ 


\— 


1~ 





1 





1++ 


1~ 


1~~ 


2 


2 


5/24 


1+- 


0"+ 


1— 





1 


1/36 


1+- 


0"+ 


1~ 


2 


1 


5/36 


1+- 


1~ 


1~~ 





1 


1/36 


1+- 


1~ 


1~- 


2 


1 


5/36 
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